Density dependent exchange contribution to dfi/dn and compressibility in graphene 
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We calculate dfi/dn (where fi = chemical potential and n — electron density), which is associated 
with the compressibility, in graphene as a function of n, within the Hartree-Fock approximation. 
The exchange-driven Dirac-point logarithmic singularity in the quasiparticle velocity of intrinsic 
graphene disappears in the extrinsic case. The calculated renormalized dfi/dn in extrinsic graphene 

on SiC>2 has the same n~z density dependence but is 20% larger than the inverse bare density of 
states, a relatively weak effect compared to the corresponding parabolic-band case. We predict that 
the renormalization effect can be enhanced to about 50% by changing the graphene substrate. 
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The band structure of graphene (a single layer of car- 
bon atoms) , by dint of its honeycomb lattice, has linear 
dispersions near the IC and K! points ("Dirac points") 
of the Brillouin zone. Recent developments in techniques 
for fabricating conducting graphene layers have thus pro- 
vided the physics community with a unique opportunity 
to study an interacting two-dimensional (2D) massless 
Dirac fermion system using table-top experimental equip- 
ment. This has led to a veritable explosion of both ex- 
perimental and theoretical activity in this field 

Around the Dirac points (which we take to be the zero 
of energy), the kinetic energy for a "bare" electron (see 
below) is = si>o|fc|, where k is the wavevector with 
respect to the Dirac point, and s = +1 and —1 for the 
conduction and valence bands, respectively. The elec- 
tron chemical potential ft, which in intrinsic graphene 
is at zero, can be shifted up or down by doping and/or 
application of external gate voltages, with a concommi- 
tant change in the electron density. This paper reports 
the calculation of dfi/dn, which is related to the elec- 
tronic compressibility, in extrinsic graphene at tempera- 
ture T = as a function of the density n. [In this pa- 
per, unless otherwise indicated, partial derivatives are at 
constant area and T = 0, and n refers exclusively to the 
free carrier density (i.e., the difference in electron density 
from that of intrinsic graphene) in the gated graphene, 



which we take to be substantially less 
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than the intrinsic electron density n v (> 10 15 cm -2 ) fill- 
ing up the valence band.] We obtain fi(n) by evaluating 
the electron self-energy within the Hartree-Fock approx- 
imation (HFA). The HFA is a good approximation up 
to reasonably high values of r s (~ the ratio of the aver- 
age carrier potential to kinetic energy) in parabolic-band 
semiconductors, and we expect it to also give reliable re- 
sults in graphene, where r s < 1. 

It is useful (and conceptually meaningful) to divide 
2D graphene into three different systems depending on 
the band filling: bare or empty, a theoretical abstraction 
of just one electron in the graphene honeycomb lattice 



as appropriate for the single-particle band-structure cal- 
culation with both valence and conduction bands com- 
pletely (and unphysically) empty, or equivalently, the un- 
physical situation where the interaction between the elec- 
trons is turned off; intrinsic, i.e., the undoped and un- 
gated situation, which is a zero-gap semiconductor with 
a completely full (empty) valence (conduction) band and 
chemical potential fi (— Ep, since we are at T = 0) 
precisely at the Dirac point; extrinsic, i.e., gated/doped 
graphene with a tunable 2D free carrier density n of elec- 
trons (holes) in the conduction (valence) band, with \i be- 
ing above (below) zero, i.e. in the conduction (valence) 
band. Note that only the empty system can be char- 
acterized by the bare, noninteracting parameters (e.g., 
velocity vq, density of states Do) with both intrinsic and 
extrinsic graphene being characterized by renormalized 
parameters. We emphasize that the bare graphene pa- 
rameters, being unphysical abstractions, cannot be ex- 
perimentally determined. 

In the absence of interaction, dfi/dn is just the in- 
verse of the bare or non-interacting single-particle den- 
sity of states at the Fermi level: (dn/dfx)o = D (Ep) = 
\/9 n ~/(\/Trhvo), where vq is the bare graphene carrier ve- 
locity associated with the linear energy dispersion, and 
g(= 4) is the product of the spin (g s — 2) and valley 
(g v = 2) degeneracies of the graphene carriers. This im- 
plies that d/i/dn in graphene is a direct measurement of 
the thermodynamic Fermi velocity renormalization due 
to electron-electron interaction effects. (This should be 
distinguished from the quasiparticle Fermi velocity renor- 
malization, as discussed later.) Our goal here is to the- 
oretically calculate the renormalized dfi/dn in intrinsic 
and extrinsic graphene including exchange interaction ef- 
fects, or equivalently in the HFA, which should be an ex- 
cellent quantitative approximation in 2D graphene. Our 
calculated carrier density dependence of dfi/dn can be 
directly compared to experimental measurements in ex- 
trinsic graphene. 



2 



The exchange self energy is given by [2] 

Zx*(k) = -J2nF(Sk- q s)Vc(q)F.*(k,k-q), (1) 

s'q 

where s,s' = ±1 are the band indicies, and since we 
assume T = 0, the fermi function n F (^) — or 1 for 
£ = e — fj, less than or greater than 0, respectively. 
V c (q) = 2ne 2 / (nq) is the bare coulomb potential (k is the 
background dielectric constant in the graphene layer). 
F ss '(k,k') = (1 + ,ss'cos#fcfc')/2 arises from the wave- 
function overlap factor, where 9kk' is the angle between 
k and k' . We assume that the valence band is cut off at 
the wavevector k c with respect to the Dirac point. The 
ultra-violet regularization associated with the wavevector 
cut-off at k c happens at a very large wavevector, compa- 
rable to the lattice wavevector; k c ~ lA . Within the 
HFA, fi = e ( °l s + £ X:S (fc F ), where k F = (4Tr\n\/g)* is 
the Fermi wavevector. 

We separate the exchange self-energy into contribu- 
tions from the intrinsic electrons, £ x nt , and the extrinsic 
carriers, £° xt . That is, £ X , S (A) = E£* (fc)+S° x *(fc), where 

K%k) = -J2Vc(q)F s ^(k,k-q): (2a) 
ZTs(k) = -J2 Sn Fi^- q , s ')V c (q)F ssl (k,k-qpb) 

s'q 

where 5n F (£ k - q y) = n F (^k-q,s') — ^0- — s') is the differ- 
ence in the electron occupation from the intrinsic T = 
case. Evaluating the integrals, we obtain 
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where =p in Eq. (|3b[) is for \i ^ 




if x < 1; 
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x[K(y)-E(y)-fy 2 ] 
xff'dy [K(y)-E(y)], 



for a; < 1; 
for x > 1. 



(5) 

Here, -K"(x) and -E(x) are the complete elliptic integral of 
the first and second kinds, respectively Note that the 
T = exchange self-energy for a regular parabolic-band 
two-dimensional electron gas (2DEG) is [H Y^°(k) = 

-Hj^fCfe)- At T = °' Ex ( fc ) does not de P end on the 
band-structure away from the Fermi surface (since n F is 
either 1 or for £ < and £ > respectively, independent 



of the details of the band-structure) and therefore the 
only difference between E x (fc) for the parabolic-band case 
and and the intraband contribution for graphene is the 
difference in the wavefunction overlap factor F ss i(k,k'). 
This accounts for presence of the f(k/k F ) in both the 
T = expressions for EP b (fc) and E x xt (A)- 

First, we examine the intrinsic self-energy, E™ (k), 
which is independent of the carrier density n. Since 
we are interested in the states around the Dirac point, 
the argument of the functions / and h in Eq. (|3a|). 
k/k c < 1. For small x, f(x) = E(x) «| - 0{x 2 ), and 
h{x) «x(f log(|)-f ) + 0(a; 3 ) (in Eq. ©, the integrand 

y as y — > 0, and therefore the 
Therefore, 



y- 3 {K(y)-E(y)-iyV 



integral 



for small 



(6) 

The term — e 2 k c /(2n) in Eq. ([5]) simply shifts energy zero 
and can be ignored. The other terms renormalize the 
quasiparticle velocity. Ignoring terms of order k/k c the 
renormalized quasiparticle velocity is 0, Q 
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where fe c = 4e~Sfc c w 0.9fc c and r s u; = e 2 / (hnvo). Exper- 
imental measurements of the quasiparticle velocity in in- 
trinsic graphene will yield Wj n t (in the absence of phonon 
coupling), and not the bare velocity vq, which applies 
only for the unrealistic situation of a completely empty 
valence-band. This situation is analogous to the quan- 
tum electrodynamics calculation of the self-energy of a 
bare electron. The bare electron charge and mass of the 
theory are never observed. Instead, experimentally one 
sees the scale-dependent renormalized charge and mass, 
which include effects of the electron self-energy. The log- 
arithmic dependence of the intrinsic graphene velocity is 
probably difficult to observe because of the smallness of 
the prefactor r s °' ) /4 w 0.2 for graphene mounted on a 
SiC>2 substrate with one side exposed to air (hence, the 
effective k in the graphene layer is the average of the k 
of air and SiC>2, ~ 2.5). To see clearly the logarithmic 
dependence in Eq. (0, k must be varied over a fairly 
wide range. Furthermore, the logarithmic divergence in 
Wi n t(A) at k — > occurs only in the intrinsic graphene, 
and not in the extrinsic case. 

In extrinsic graphene, k F ^ 0. For k/kp <C 1, the 
small x expansions for f(x) and h(x) in Eq. (|3b[) yield 
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(where =F is for fi ^ 0). The log(A) term in Y££l{k — > 0) 
cancels the equivalent term in Yj™ A s (k — > 0), so the sum, 
S XjS (A) 



£ x xt (A) + ££*(Jfe), has a finite derivative at 
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FIG. 1: (Color online) Exchange self-energies for graphene 
(with (i > 0), and for a parabolic-band 2DEG, as functions of 
wave vector. Note that £ x +(fc) + E^(fc) = £ x b (fc)- 



k = 0, and the renormalized velocity in the extrinsic 
1 



, has no fc — * 



case, w GX t(fc = 0) = v a 
logarithmic divergence. 

Fig. 1 shows £ x xt (fc) for graphene with /i > and, for 
comparison, £P b (fc). At k — kp (and fj, > 0), 



K^kp) = - 



1 



"-5 



Sf(fcF), (9) 



where C w 0.916 is Catalan's constant, and (fcjr) = 
— ire 2 kp / (ttk). As in the case of £P b (fc), the slope of 
^x+(^) f° r graphene with fi > has a logarithmic di- 
vergence as k — > fcf from both the /(x) and terms 
in Eq. (|3b|) . We expect that the logarithmic divergence 
in dT^/dk will disappear when correlation effects are 
included, as in the case of the parabolic-band £P b . Note 
that this logarithmic divergence has no singular patho- 
logical effect on dfj,/dn, the quantity of interest in this 
work, and is irrelevant for our purpose. The has a 
finite derivative at k — kp, because for an electron-doped 
sample there is no Fermi surface at k = kp in the valence 
band. 

We now discuss the effect of the exchange self-energy 
on dfj,/dn. In a regular parabolic-band 2DEG with mass 



in. 



\dnJ Q 



2h 2 n 



is density independent. When inter- 

actions are included this picture changes due to exchange 
and correlation effects of the Coulomb potential. Within 

the HFA, which neglects correlation effects, ff^J = 

(ilj b ( X _ ^^ b ), where = V2e 2 m/{nk F ) oc n~*. 
Thus, d/i/dn becomes negative at low enough densities. 
Measurements of d/i/dn in two-dimensional electron and 
hole gases have confirmed this behavior 0, The ob- 
served change of sign in d[ijdn comes mainly from the ex- 
change contribution to the total energy. It is known that 
the correlation corrections to djji/dn beyond the HFA is 



not very large (< 20%) [§], even at the reasonably large 
effective r s (3> 1) values at which the 2D semiconductor 
experiments have typically been carried out. 

What is the contribution of exchange on d/i/dn in 
graphene? Using kp = y^r \n\, together with (j, s (kp) — 
M j_ E^Kkr) + ^Ts(k F ), and Eq. © gives 



k F ,s 
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where 



e /(kHv™ 1 ) [here, v 



int 



(10) 



t(fcp)] and 



(d[i/dn)\ n t = (ftw^ lt v / 7r)/(2 1 /|n|) is the inverse of the 
density of states for intrinsic graphene. This shows 
that in extrinsic graphene, the exchange effect changes 
dfi/dn from the bare value (d[i/dn)o by a factor of 



,(°) 



>r from the intrinsic 
value {dfi/dn)i n t by a factor of — 7r _1 (C + \)r™ 1 . Thus, 

(9/z/9n) oxt is enhanced over (d/j/dn)o by « (0.25) ri°^ 
(for n = 10 12 cm~ 2 ), but reduced from (d/j,/ dn)i n t by 
a factor of f=a (— 0.45) r^. 11 *. If we take single-particle 
band-structure graphene velocity vq w 10 8 cm/s, we get 
0.9 for Si02 mounted graphene, giving an ex- 



change enhancement of approximately 20% over the bare 
(d/j,/dn)o; see Fig. 2(a). Estimating the change with re- 
spect to (d/j,/dn)i n t is a little trickier because r™*, which 
depends on the intrinsic graphene velocity Vp*, is un- 
known since the intrinsic graphene velocity is at present 
unknown! An approximate way to estimate the instrin- 
sic r* nt is to change k — > kk* where k* is the effect of 
the background screening by the filled valence band[l0(, 
with k* = 1 + Zg s g v ri 0) « 2, which gives rf* w rf ] /2, 
leading to around a 20% decrease of (d^/dn) cxt with re- 
spect to ((9/x/9n)int- The extrinsic dfi/dn depends on k, 
as shown in Fig. 2(b), but unlike the parabolic band case, 
for experimentally relevant parameters (d^i/dn) oxt does 
not change sign. By using freely suspended graphene 
(i.e., k — 1), the many body renormalization can be en- 
hanced to around 50%. 

In the literature, dfi/dn is often associated 
with the term "compressibility," defined as 
K = —V~ 1 (dV/dP)T,N, where N is particle num- 
ber, V is the system volume/area, P = ~{dT / 8V)t.,n is 
the pressure, and T is the Helmholtz free energy. It can 
be shown that 11] K^ 1 = n 2 (dfi/dn), where n = N/V. 
In experimental papers on the compressibility of electron 
gases, the quantity that is measured is not actually 
the compressibility (after all, experimentalists do not 
physically compress the electron gas and measure the 
change in pressure) but dfi/dn, which is then converted 
to K^ 1 by multiplication of n 2 . In graphene, it is in 
fact ambiguous which n should be used — (a) the free 
carrier density or (b) the density of the electrons in the 
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FIG. 2: (Color online) (a) Calculated dfi/dn as a func- 
tion of free carrier density, using the following parameters: 
k c = 1/a (a =2.46 A), k = 2.5, v = 10 8 cm/s, and ri 0) = 0.9. 
The "bare" curve is dfi/dn of a noninteracting graphene, and 
the "int" and "ext" curves are for the intrinsic and extrinsic 
cases, respectively. The inset shows the ratios of the renor- 
malized intrinsic velocity at &f and the inverse of the extrin- 
sic compressibility to their corresponding bare quantities, (b) 
(d[i/dn) e xt for different values of k (hence, different r s ). 



band (i.e., free carrier density plus n v )7 The answer 
depends on which hypothetical compressibility is being 
considered — (a) corresponds to one in which the area 
enclosing the free carriers is changed but the underlying 
lattice is kept constant, and (b) to one where the volume 
of the underlying lattice (and hence n v ) also changes. 
To avoid any ambiguities, we use the quantity dfi/dn. 

Before concluding, we point out that it is incorrect 
to think of (<9/x/<9n) cxt as providing a measurement of 
many-body quasiparticle Fermi velocity renormalization 
in graphene. In particular, the quasiparticle velocity 
renormalization is given by the Ff parameter in Fermi 
liquid theory through Galelian invariance vp/v r p n = 
1 + F*, whereas the renormalization of dfi/dn is related 
to the Fermi liquid parameter Fq through the identity 
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dn 
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Hence, although (dfi/dnjo for the bare system is propor- 
tional to the bare particle velocity vq at the Fermi sur- 
face, (dn/dn) ext for the extrinsic case is not proportional 
to the quasiparticle velocity because of the presence of 
the additional Fermi liquid parameter Fq. 

We conclude with a discussion of the possible effects 
of disorder and correlation on graphene dnjdn. Wc 
believe that correlation, neglected in our Hartree-Fock 
theory, would introduce only small quantitative correc- 
tions to our calculated results, particularly because of 
the relatively small values of r s (< 1) in graphene. 
Thus, our Fock exchange approximation for graphene 
dnjdn should quantitatively be an excellent approxima- 
tion. Disorder would also introduce only small quanti- 
tative corrections except at low extrinsic carrier densi- 
ties (|n| <5x 10 11 cm -2 ) associated with the so-called 
"minimal graphene conductivity" regime, where random 
charged impurities in the substrate introduce [12|, |13| in- 
homogeneous electron/hole puddles in the graphene layer 
which would lead to random spatial variations in dfi/dn 
over 5-20 nm (10 - 100 meV) length (energy) scales. 
Finally, finite temperature would have little effect on our 
results because Ep > 1000 K in the usual density range 
of experimental interest (>5x 10 11 cm -2 ). 

This work is supported by US-ONR and LPS-NSA. 

Note added: After submission of this manuscript, we 
received a preprint [3] reporting an experimental obser- 
vation of the exchange contribution to dnjdn that is con- 
sistent with our theory. 
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